UNCLASSIFIED 


_ AD  NUMBER _ 

AD112418 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited.  Document  partially  illegible. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
and  their  contractors ; 

Administrative/Operational  Use;  07  DEC  1956. 
Other  requests  shall  be  referred  to  Deputy 
Chief  of  Staff  Research  and  Development , 
Department  of  the  Air  Force,  Washington,  DC 
20330.  Document  partially  illegible. 


_ AUTHORITY 

Rand  per  ltr,  20  Jun  1969 


THIS  PAGE  IS  UNCLASSIFIED 


firmed  Services  Technical  Information  Agency 

Reproduced  by 

DOCUMENT  SERVICE  CENTER 

KNOTT  BUILDING,  DAYTON,  2,  OHIO 

This  document  is  the  property  of  the  United  States  Government.  It  is  furnished  for  the  du¬ 
ration  of  the  contract  and  shall  be  returned  when  no  longer  required,  or  upon  recall  by  ASTIA 
to  the  following  address:  Armed  Services  Technical  Information  Agency, 
Document  Service  Center,  Knott  Building,  Dayton  2,  Ohio. 


NOTICE:  WHEN  GOVERNMENT  OR  OTHER  DRAWINGS,  SPECIFICATIONS  OR  OTHER  DATA 
ARE  USED  FOR  ANY  PURPOSE  OTHER  THAN  IN  CONNECTION  WITH  A  DEFINITELY  RELATED 
GOVERNMENT  PROCUREMENT  OPERATION,  THE  U.  S.  GOVERNMENT  THEREBY  INCURS 
NO  RESPONSIBILITY,  NOR  ANY  OBLIGATION  WHATSOEVER;  AND  THE  FACT  THAT  THE 
GOVERNMENT  MAY  HAVE  FORMULATED,  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE 
SAID  DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
IMPLICATION  OR  OTHERWISE  AS  IN  ANY  MANNER  LICENSING  THE  HOLDER  OR  ANY  t-TPER 
PERSON  OR  CORPORATION,  OR  CONVEYING  ANY  RIGHTS  OR  PERMISSION  TO  MANU.  CtURE, 
USE  OR  SELL  ANY  PATENTED  INVENTION  THAT  MAY  IN  ANY  WAY  BE  RELATED  THERETO. 


ASS 


Best 

Available 

Copy 


RESEARCH  MEMORANDUM 


This  is  a  working  paper.  It  may  be  expanded, 
modified,  or  withdrawn  at  any  time. 
The  views,  conclusions,  and  recommendations 
expressed  herein  do  not  necessarily  reflect  the 
official  views  or  policies  of  the  U.  S.  Air  Force. 


k-fl  FI  (fafmttfiott 


SANTA  MONICA 


CALIFORNIA 


U  S  AIR  FORCE 

PROJECT  RAND 

RESEARCH  MEMORANDUM 


NOTES  ON  LINEAR  PROGRAMMING:  PART  XXXVI: 
THE  ALLOCATION  OF  AIRCRAFT  TO  ROUTES  — 
AN  EXAMPLE  OF  LINEAR  PROGRAMMING  UNDER 
UNCERTAIN  DEMAND 


Allen  R.  Ferguson 
George  B.  Dantzig 


RM-1833 


ASTIA  Document  Number  AD  112418 
7  December  1956 


Auigned  lo 


This  is  a  working  paper.  It  may  be  expanded,  modified,  or  with¬ 
drawn  at  any  time.  The  views,  conclusions,  and  recommendations 
expressed  herein  do  not  necessarily  reflect  the  official  views  or 
policies  of  the  United  Slates  Air  Force. 


nit 


WMD 


WOO  MAIN  IT. 


SANTA  MONICA 


CAllFOtNIA 


Copyright  1957 
The  RAND  Corporation 


RM-1833 
1 2-7-5  o 
-li- 


SUMMARY 

The  purpose  of  this  paper  is  to  illustrate  an  application 
of  linear  programming  to  the  problem  of  allocation  of  aircraft 
to  routes  in  order  to  maximize  expected  profits  when  there  is 
uncertain  customer  demand.  The  approach  is  intuitive;  the 
theoretical  basis  of  this  work  is  found  in  an  earlier  study. 

The  allocations  are  compared  with  those  obtained  under  the 
usual  procedure  of  assuming  a  fixed  demand  equal  to  the  expected 
value.  The  computational  procedure  is  similar  to  that  of  the 
fixed-demand  case,  with  only  slightly  more  computational  effort 
required . 

This  paper  is  intended  both  for  readers  interested  in 


routing  problems  (and  analogous  resource-allocation  problems) 
and  for  those  interested  in  studying  an  example  of  an  appli¬ 
cation  of  linear  programming  under  uncertainty. 
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THE  ALLOCATION  OF  AIRCRAFT  TO  ROUTES— AN  EXAMPLE 
OF  LINEAR  PROGRAMMING  UNDER  UNCERTAIN  DEMAND 

1  .  INTRODUCTION 

There  are  many  business,  pnonomic,  and  military  problems 
that  have  the  following  characteristics  in  common:  a  limited 
luantity  of  capital  equipment  or  final  product  must  be  al¬ 
located  among  a  number  of  final-use  activities,  where  the  level 
of  demand  for  each  of  these  activities,  and  hence  the  payoff,  is 
uncertain;  further,  once  the  allocation  is  made,  it  is  not  eco¬ 
nomically  feasible  to  reallocate  because  of  geographical  sepa¬ 
ration  of  the  activities,  because  of  differences  in  form  of  the 
final  products,  or  because  of  a  minimum  lead  time  between  the 
decision  and  its  Implementation.  Examples  of  3uch  problems  are 
(l)  the  scheduling  of  transport  vehicles  over  a  number  of  routes 
to  meet  a  demand  in  some  future  period  and  (2)  the  allocating  of 
quantities  of  a  commodity  at  discrete  time  intervals  among 
several  storage  or  distribution  points  while  the  future  demand 
for  the  commodity  is  unknown.  It  is  assumed,  however,  that 
demand  can  be  forecast  or  estimated  as  a  distribution  of  values, 
each  with  a  specified  probability  of  being  the  actual  value. 

The  general  area  where  the  techniques  of  this  paper  apply 
may  be  schematized  broadly  as  problems  where: 

1.  Alternative  sets  of  activity  levels  can  be  chosen 
consistent  with  given  resources. 

2.  Each  set  of  chosen  activity  levels  provides  the 
facilities  or  stocks  to  meet  a  demand  that  is  itself 
unknown  but  that  has  a  known  frequency  distribution. 
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') .  Profits  depend  on  the  costs  ot  the  facilities,  or 
stocKs,  ami  on  the  revenues  from  the  demand. 

a.  The  reneral  objective  is  to  determine  that  set  of 
active  sy  level  a  that  mav  1  -• !  /.os  profits. 

The  pan'T  ‘nil  Lied  "Linear  Programming  under  Uncertainty" 
[,J  \-j rn.s  tne  theoretical  basis  for  the  present  work.  Our  pur¬ 
er;,;'  is  to  i.  1 U  tra t<’  tne  (  rocedura  l  steps  wltn  the  example 
tint,  in  i  set,  rl,;  Iris  l  ly  prompted  the  referenced  theoretical 
work  1  ri  this  uva.  Thus  ,  1 1 1 1 1  < •  Ln  tne  way  of  rirorous  theory 
will  be  attempt...1 1  In  tills  papier,  althoupn  each  step  will  be 
.instil  l':  ;  Intuit  ively  . 

Tne  method  Is  explained  by  the  use  of  a  mode  j.  for  routinr 
aircraft.  Tevral  types  of  aircraft  are  allocated  over  a  number 
of  rout  os ;  s :  monthly  .lemand  for  servi  e  over  esen  route  is 
assumed  to  us  mown  on  Ly  as  a  d  Is L rl  butl.on  of  probable  vs  l.u  .  . 

Tr.>  a  .1.  r\;  ra  f  t.  <•-.»* :  r»o  a  1  loci  led  an  to  minimize  the  sun.  of  (a)  the 
ei.'S t  of  perform iu  r,  th  :  transportation  and  (b)  tne  ex|  e  te  1  value 
of  tne  reveiiu 1  lost  through  the  failure  to  serve1  all  the  Lraffl 
t;  at  us  lus  'l  .1 y  u  :  vs:  j  oj  s  . 

For  .  in';  of  month- -to  month  s  :rie  lu.L  I  n /  ,  m  '  1  r  tr  insport 

o:  •  r  1  tor  /.ot.;  ..  p  r  ■■su.miuly  ,  1  ■  1 1  •  ■  r‘  iOouL  usvin  ,  tj  :.  .  ■  w. 
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It  is  suggested  that  the  reader  make  sensitivity  tests  by  modi¬ 
fying  the  demand  distributions  given  in  the  illustrative  example. 

Passenger  demand,  of  course,  occurs  on  a  day-by-day — in 
fact,  on  a  fllght-by-f light  basis.  The  assumed  number  of  pas¬ 
sengers  per  aircraft  of  a  given  type  per  flight  on  a  given  route 
may  be  thought  of  as  an  ideal  number  that  can  be  increased 
slightly  by  decreasing  the  amount  of  air  freight  when  this  is 
Indicated,  and  by  "smoothing"  the  demand  through  encouraging 
the  customers  to  take  open  reservations  on  alternative  flights 
as  opposed  to  less  certain  reservations  on  desired  flights.  In 
spite  of  these  possible  adjustments,  traveler  preferences  and 
the  inevitable  la3t-minutc  cancellations  do  cause  loss  of  pas¬ 
senger-carrying  capability.  However,  the  best  way  to  reflect 
these  effects  of  the  daily  variations  in  demand  are  beyond  the 
3cope  of  this  paper.  For  our  purpose  here,  either  the  aircraft 
passenger-carrying  capability  or  the  demand  may  be  thought  of  as 
adjusted  downward  to  reflect  the  loss  due  to  daily  variations 
of  demand . 

The  method  employed  is  simple,  and  the  example  used  can  be 
solved  by  hand  in  an  hour  or  two.  Larger  problems  can  be  solved 
with  computing  machines. 

In  a  previously  published  paper  [l]  ,  the  method  was  applied 
to  the  same  example,  assuming  the  demand  on  each  route  to  be 

known;*  the  present  paper  continues  the  analysis,  showing  how 

- 

This  was  equivalent  to  using  the  expected  value  of  demand, 
rather  than  taking  account  of  the  whole  frequency  distribution, 
as  in  the  present  paper. 
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1.0  handle  a  frequency  distribution  of  demand  relative  to  each 
route.  A  different  allocation  is  found  to  be  optimal  in  this 

case  . 

We  shall  row  describe  the  problem,  briefly  indicate  the 
nature  of  the  solution  based  on  expected  values,  show  the  method 
of  solving  the  problem  using  stochastic  values  for  demand,  and 
finally  compare  the  two  solutions. 
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2.  REVIEW  OF  FI  XED-DEMAND  EXAMPLE 

The  fixed-demand  example  that  we  are  using  to  illustrate 
the  method  takes  a  fixed  fleet  of  four  types  of  aircraft,  as 
shown  in  Table  1 . 


Table  1 

ASSUMED  AIRCRAFT  FLEET 


Type 

Description 

Number 

Available 

A 

Postwar  4-englrie 

10 

B 

Postwar  2-engine 

19 

C 

Prewar  2— engine 

25 

D 

Prewar  4-engine 

15 

These  aircraft  have  differences  in  speed,  range,  payload  capacity, 
and  cost  characteristics.  The  assumed  routes  and  expected  traffic 
loads  (the  distribution  of  demand  will  be  discussed  later)  are 
given  in  Table  2. 
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Table  2 


TRAFFIC  LOAD  BY  ROUTE 


Route 

Route 

Miles 

.  ..  . 

Expected 
Number  of  , 
Passengers0 

Price 
One-way 
Ticket  ($) 

(1) 

N.Y.-L.A . 

( 1-stop ) 

2,473 

o 

o 

o 

JT  \ 
CM 

130 

(2) 

N.Y.-L.A . 

(2-stop) 

2,47b 

12,000 

130 

(3) 

N . Y .—Dallas 

(0-stop ) 

1,381 

18,000 

70 

(4) 

N . Y .-Dallas 

( 1-stop ) 

1 ,4  39 

9,000 

70 

(o) 

N  .Y  .-Boston 

( O-o top ) 

1 8b 

00, 000 

10 

'lOH'icial  Airline  Guide,  July,  1934,  p.  270.  The  New  York- 
Los  Angeles  routes  are  via  Chicago  and  via  Chicago  and  Denver; 
the  atop  on  route  between  New  York  and  Dalian  ia  at  Memphi3 . 

°Thls  la  the  expected  number  of  full  one-way  trip3  per 
month  to  be  carried  on  each  route.  If  a  passenger  get3  off 
en  route  and  la  replaced  by  another  passenger,  it  is  counted 
as  one  full  trip. 


Since  this  paper  proposes  to  Illustrate  the  applicability 
of  a  method  of  solving  problems  in  which  several  realistic  ele¬ 
ments  are  considered,  it  is  assumed  that  not  all  aircraft  can 
carry  their  full  loads  on  all  routes  and  that  the  obtainable 
utilisation  varies  from  route  to  route.  Specifically,  Type  B 
Is  assumed  to  be  able  to  operate  at  only  7b  per  cent  payload  on 
Route  3,  and  Type  D  at  80  per  cent  on  Route  1;  Type  C  cannot  . 
fly  either  Route  1  or  Route  /),  and  Type  B  cannot  fly  Route  1. 
Utilization  is  defined  as  the  average  number  of  hours  of  useful 
work  performed  per  month  by  each  aircraft  assigned  to  a  particular 
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route .  Utilization  of  500  hours  per  month  is  assumed  on  Routes 
1  and  2,  285  on  Routes  5  and  4,  and  240  on  Route  5* 

The  assumed  dollar  costs  per  100  passenger-miles  are  shown 
in  Table  5-  These  do  not  include  any  capital  costs  such  as  those 
of  the  aircraft  and  ground  facilities.  They  represent  variable 
costs  such  as  the  cost  of  gasoline,  salaries  of  the  crew,  and 
costs  of  servicing  the  aircraft. 

A  second  sort  of  "cost"  is  the  loss  of  revenue  when  not 
enough  aircraft  are  assigned  to  the  route  to  meet  the  passenger 
demand.  In  this  case,  the  loss  of  revenue  is  the  same  as  the 
price  of  a  one-way  ticket  shown  in  the  E  row  of  Table  5. 


Table  5 

DOLLAR  COSTS 


Type  of 
Aircraft 

Route 

1  -  N.Y. 
to  L.A . 
l-stop  ($) 

2  -  N.Y. 
to  L.A. 
2-stop  ($) 

5  -  N.Y. 
to  Dallas 
0-stop  ($) 

4  -  N.Y. 
to  Dallas 
1-stop  ($) 

5  -  N.Y. 
to  Boston 
0-stop  ($) 

Per  100  Passenger-miles 

1  -  A 

0.45 

0.57 

0.45 

0.47 

0.64 

2  -  B 

— 

0.54 

0.85 

0.65 

0.88 

5  -  C 

— 

0.92 

— 

0.95 

1.15 

4  -  D 

0.74 

0.51 

0.59 

0.52 

0.8l 

Per 

3. 

Passenger  Turned  Away 

5  -  E 

1 

150 

150 

70 

70 

10 

(u) 

(O) 

(7) 

(7) 

(1) 

Figures  shown  in  parentheses  are  1000' s  of  dollars  lost  per 
100  passengers  turned  away.  (Throughout  this  paper,  passengers 
are  measured  in  units  of  hundreds.) 


Baaed  on  the  speeds,  ranges,  payload  capacities,  and  turn¬ 
around  times,  passenger-carrying  capabilities  were  determined. 

The  resultant  potential  number  p^j  (in  hundreds)  of  passengers 
that  can  be  flown  per  month  per  aircraft  of  type  i  on  Route  J 
is  shown  in  Table  4;  see  the  staggered  upper  right  figure  in 
each  box.  By  multiplying  these  numbers  by  the  corresponding 
costs  per  TOO  passenger-miles  given  in  Table  d  and  by  the  number 
of  miles  given  in  Table  2,  the  monthly  coot  per  aircraft  can 
also  be  obtained.  This  is  given  in  the  lower  left  figure  c^  , 

Ln  each  box;  explicitly,  c^.  is  the  cost  in  thousands  of  dollars 
per  month  per’  aircraft  of  type  i  assigned  to  the  Route  J.  The 
re' venue  losses  e,  jf  In  thousands  of  dollars  per  100  passengers 

not  carried,  are  given  ln  the  E  row  of  Tabic  4;  finally,  we  do- 
* 

fine  p(  ,  =  1 .  The  staggered  layout  of  Table  4  was  chosen  so 

J  d 

as  to  identify  the  correspond  In,'  data  found  in  Table  p;  the 
latter  is  the  work  sheet  upon  which  the  entire  problem  is  solved. 

The  basic  problem  is  that  of  determining  the  number’  of  air¬ 
craft  of  each  type  to  assign  to  each  route  consistent  with  air¬ 
craft  availabilities  (Table  1 )  and  of  determining  how  much  rev¬ 
enue  will  be  lost  due  to  failure  of  allocated  aircraft  to  meet 
passenger  demand  on  various  routes  (Tables  2  and  j>)  .  Since  many 
alternative  allocations  are  possible,  our  specific  objective  will 
be  to  find  that  allocation  that  minimizes  total  costs,  where 
costs  are  defined  as  operating  costs  plus  lost  revenues  based 

on  the  cost  factors  given  in  Table  ). 

~  . 

This  will  make  it  easier- to  form  the  passenger-balanc  -  or 
"column"  equations  (c1). 
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Table  4 

PASSENGER-CARRYING  CAPAPILITIESa  AND  COSTSb 


Route 

Type  of 
Aircraft 

1  -  N.Y. 
to  L.A. 

1 -stop 

2  -  N.Y  . 
to  L.A. 
2-stop 

a  -  N.Y.  | 
to  Dallas 
O-stop 

4  -  N.Y. 
to  Dallas 
1-stop 

5  -  N.Y 
to  Boston 
0-stop 

Per  Aircraft  per  Month 


1  -  A 

Pi i=16 

c  1 1  =  3  8 

Pi  2=  16 

c  1 2  =  2  1 

P  1  4=25 

C  1  4  =  10 

Pis=8i 

c  1  5=  10 

2  -  B 

-# 

P  2  2  =  1 0 

P23=l4 

P  2  4  =  1 5 

P25=57 

c22=16 

c  2  3  =  1 6 

c  2  4  =  1 4 

0  2  5  =9 

5  -  C 

* 

P  3  2=5 

* 

P34=  ( 

P35=29 

e  3  2= 1 0 

c  34  =  9 

C  35  =  0 

4  -  D 

P4  1=9 

P-V  2=  1  1 

P43=22 

1 

P  4  4  =  1  7 

mm 

c4 i=17 

c42=16 

c43  =  17 

0  4  4  =  15 

Per  100 

Passengers 

Not  Carried 

(Losses ) 

5  -  E 

Psi  =  l 

P  5  2  =  1 

P53=l 

P54  =  l 

P55  =  l 

C  5  1  =  13 

*-52  =  15 

C53=  I 

*-5  4  =  7 

055  =  1 

a 

b 


Capabilities  p^j  are  measured  in  hundreds  of  passengers. 
Costs  are  measured  in  thousands  of  dollars. 
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Th  1 3  fixed  -demand  model  may  be  formulated  inathemat  Lcaily 
as  a  linear  propramrnln;;  problem.  Let  x^  (  denote  the  unknown 
number  of  aircraft  of  tne  1^  type  assigned  to  J^1  route,  where 


n-1 .  If  x.  denotes  the 
’  in 


th 
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FIXED-DEMAND  MODEL 
Find  numbers  x^(,  and  the  minimum  value  of  z,  satisfying 
the  following  conditions. 

(1)  How  Guns:  x^  +  +  •••  +  x^n  =  a^ 

( J  -  1 ,  2  ,  .  .  . ,  m ) , 

(2)  Column  Sums :  [>]y<]  ,  t-  +  +  P,njxmj  =  dj 


Any  set  of  assignments  x^  ,  satisfying  Eqs .  (l),  (2),  and  (5) 
is  termed  a  feasible  solution ,  and  a  feasible  choice  that  min¬ 
imizes  the  total  cost  z  of  tne  assignment  is  called  an  optimal 
( feasible  )  solution  . 

Table  5  shows  the  optimal  assignment  of  aircraft  to  routes, 
based  on  fixed  demand,  a3  developed  in  the  earlier  study.  The 
values  assigned  to  the  unknowns  x^j  appear  underlined  in  the 
upper  left  of  each  box  unless  x^j  =  0  in  which  case  it  is  omitted; 
the  entire  layout  takes  the  form: 


c 
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The  suras  by  rows  of  the  x^  ^  entries  in  Table  b  equated  to  avail¬ 
abilities  yield  Eqs  .  (l).  The  sum  by  columns  of  the  ,  weighted 
by  corresponding  values  of  p^  equated  to  demands  yield  Eqs .  (2); 

'•he  x,  ,  weighted  by  corresponding,  c.  and  summed  over  the  entire 
■L  .1  i  J 

'.able  yield  Eq.  (4).  As  noted  earlier,  Table  b  is  actually  the 
work  sheet  upon  which  the  entire  problem  is  solved.  Later  we 
shall  discuss  a  revision  of  this  work  sheet  for  solving  problems 
with  variable  demand.  All  figures  in  the  table,  except  for  the 
upper  left  entries  x^ j  and  values  of  the  so-called  "implicit 
prices"  u^  and  v,  shown  in  the  margins,  are  constants  that  do 
not  chang, e  during  the  course  of  computation.  The  values  of  the 
variables  x^,  u^,  and  v^,  however,  will  change  during  the  course 
of  successive  iterations  of  the  simplex  method  a3  adapted  for 
this  problem.  For  this  reason  it  is  customary  to  cover  the  work 


sheet  with  clear  acetate  and  to  enter  the  variable  information 


with  a  grease  pencil  so  that  the  marks  can  be  easily  erased; 
alternatively,  a  blackboard  or  semitransparent  tissue-paper 
overlays  can  be  used.  The  detailed  rules  for  obtaining  the 
optimal  solution  shown  are  given  in  [l"J  and  will  not  bo  repeateu 
here.  Instead,  a  more  general  set  of  rules  for  the  uncertain- 
demanu  case  will  be  given;  these,  of  course,  could  be  used  in 


particular  for  the  expected-demand  ease. 

In  the  following  outline  we  have  a  convenient  summary  that 
serves  to  identify  and  define  the  numerical  data  entered  in 
Table  b  and  to  give  the  test  for  optimality. 


Table  ') 
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OPTIMAL  ASSIGNMENT  FOR  FIXED  DEMAND 
Operating  Costa  and  Lost  Revenues  =  $1,000,000 


Route 


Type 

of 

Air¬ 

craft 

0) 

N.  Y. 

to 

L.  A. 

1-3  top 

(2) 

N.  Y. 

to 

L.  A  . 

2-stop 

(3) 

N.  Y. 
to 

Dallas 

0-stop 

CO 

H.  Y. 

to 

Dallas 

1-stop 

(8) 

N.  Y. 

to 

Boston 

0-stop 

(6) 

Sur¬ 

plus 

Air¬ 

craft 

Air¬ 

craft 

Avail¬ 

able 

Im¬ 

plicit 

Prices 

ui 

j.  . .  | 

10=ai 

(1)  A 

WM 

28 

28 

81 

c 

IllL-J 

18 

16 

10 

0 

-171  _ 

8 

8 

6 

19  =  a2 

(2)  B 

*  * 

10 

14 

18 

87 

0 

18 

16 

14 

9 

0 

-  61 

7-8 

17.2 

28=a3 

0)  C 

*  * 

8 

*  * 

7 

29 

0 

10 

.  _  9  : 

6 

0 

-  28 

91 

8 

19=a4 

CO  0 

mm 

11 

22 

17 

88 

0 

«9 

16 

_ 17 . 

15  : 

10 

0 

-  89 

\ 

mm 

¥  # 

(5)  E 

1 

1 

1 

1 

1 

0 

Deficit 

13 

13 

7 

7 

1 

0 

0 

Demand 

6 

280 

120 

180  ! 

| 

90 

600 

* 

¥ 

Im¬ 

plicit 

Prices 

1 1  .8 

0.6 

4.8 

4  .33 

1 

\ 

0 

v 


i 


SUMMARY 


hm  \tty 

i;-  Y-Ub 

i'i  - 


f'wii:-. tant.G :  a.  =  number  of  available  aLreraft  of 

type  i 

d.  -  cx.pcct.od  par. Gender  demand  l.n  100'  .; 

|  ‘  per  month  on  Route  ,! 

i 

p'j,  -  paGocn^er-carrylnp  capability  In 
100'g  per  month  per  a  A  tv.  ra  I't  of 
type  .1.  nGGif'nod  to  Route  j 
(p  .  ,  -  1  by  do  Tin  it  Lon) 

e.  .  --  e o g t g  Ln  10001.',  of  doilara  per 
'  month  per  aircraft  of  type  i 

aoGlp.ned  to  Route  j  (c  ,  ,  Lg 

,  i  nn  1111  i 

pe l  100  paGGf.Mi^erG 

turned  away) 


a.  .  iintrAea:  x,  .  ^  number  of  aircraft  of  type  i 

aGGlrned  to  Route  J  (x  ,  .  io  100 

of  paGaen^era  turned  ’  1 
away ) 

Omitted  x.  ,  bntriea  :  x.  ,  =0  if  upper  left  entry  in  box  J.g 
iJ  1,1  miGGlnr, 


Impiicit  L'rleeo  :  u^  and  v.  are  determined  Guch  that 

u,  +  p,  ,v,  =  c,  ,  for  (i,,J)  boxer,  with 
i  i.l  a  lj 

j  >  0 — i.e.,  with  underlined  entrica . 

Note :  u  -  =  v  =0 
ith-  L  n 


Teat  for  Optimality:  Solution  is  optimal  if,  for  all  (i,j), 

the  relation  u^  +  Pj_jvj  < 
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3.  EXTENSION  OF  EXAMPLE  TO  UNCERTAIN  DEMAND 

Up  to  thi3  point  the  problem  is  identical  with  that  de¬ 
scribed  and  solved  in  our  previous  paper.  Now,  to  introduce 
the  element  of  uncertain  demand,  we  assume  not  a  known  (expected) 
demand  on  each  route  but  a  known  frequency  distribution  of  demand 
The  assumed  frequency  distributions  are  shown  in  Table  6.  Thus 
on  Route  5  (N.Y.  to  L.A.  -  2-stop)  either  5,000  or  15,000  pas¬ 
sengers  will  want  transportation  during  the  month,  with  prob¬ 
abilities  50  or  70  per  cent  respectively.  The  assumed  traffic 
distributions  are,  of  course,  hypothetical  to  illustrate  our 
method.  The  demand  distributions  on  the  five  routes  vary  over 
wide  ranges  and  have  different  characteristics;  Route  1  is  flat, 
Route  2  is  U-shaped,  Routes  5,  4,  and  5  are  unlmodular  but  have 
differing  degrees  of  concentration  about  the  mode.  Route  4  ha3 
a  distribution  with  a  very  long  tail  that  may  reflect  a  real¬ 
istic  traffic  situation. 
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Table  b 


ASSUMED  DISTRIBUTION  OF  PASSENGER  DEMAND 
>\h  ,  =  Probability  of  Demand 


Route 

Passengers 
(in  hundreds) 

Approx .  Mean 
(in  nundreds) 

Probability 
of  Passenger 
Demand  • 

Probability 
of  Equaling 
or  Exceeding 
Demand 

200  =  d i i 

0.2  =  am 

i.o  rM 

220  =  d2 i 

0.03  -  A2i 

0.8  •  If 2  j 

1  -  A 

230  =  d3 1 

230 

0-3b  “  ^3i 

j 

o .  !\  -  y3, 

270  =  d41 

0.2  -  Asl  ; 

On*  -  <14! 

300  =  dsi 

0.2  =  j  | 

0.2  -  ?51 

30  =  d ,  2 

i 

0.3  =  \2  | 

1.0  ^12 

2  -  B 

ioO  =  d22 

120 

0.7  =  Aaa 

0  .  !  ~  %  2  2 

140  =  d  1 3 

0.1 

1  .0  =  0-  3 

1 bO  =  d23 

0.2  -  A 23 

0-9  =  0  2  3 

3  -  C 

1 80  -  d3 3 

180 

0.4  =  A  3  3 

0.1  --  *33 

200  =  d43 

0.2  ~  A 43 

1  03  "^43 

220  =  dS3 

0.1  -  A.  53 

1  0.  -  *53 

1 

i 

10  =  d  14 

0.2  -  A,4 

:.o  -*14 

30  =  d  2  4 

0.2  =  A24 

|  0  , 8  2  4 

4  -  D 

80  =  d34 

90 

■0.3  = 

0.0  ~  *34 

100  =  d44 

0.2  =  *44 

=  ^44 

340  =  d54 

0.1  =  *54 

0.)  = 

880  =  d  15 

0.1  =  '\  5 

1.0  =  *15 

3  -  E 

bOO  =  d2s 

D00 

0.8  =  A  2  5 

0-9  =  ^25 

b20  =  d35 

0.1  =  A35 

0.1  =*35 
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To  Illustrate  the  essential  character  of  the  linear-pro¬ 
gramming  problem  for  the  case  of  uncertain  demand,  let  us  focus 
our  attention  on  a  single  route — say,  Route  1 — with  probability 
distribution  of  demand  as  given  in  Table  6.  Let  us  suppose  that 
aircraft  assigned  to  Route  1  are  capable  of  hauling  1.00  Yj  pas¬ 
sengers.  The  first  200  units  (in  hundreds  of  passengers)  of 
this  capability  are  certain  to  be  used,  and  revenues  from  this 
source  (negative  costs)  will  be  13  =  units  (in  thousands  of 
dollars)  per  unit  of  capability.  The  next  20  units  of  thl3 
capability  will  be  used  with  probability  ^2i  =  0.8.  Indeed, 

80  per  cent  of  the  time  the  demand  will  be  220  units  or  greater, 
while  20  per  cent  of  the  time  it  will  be  200  units;  hence,  the 
expected  revenues  per  unit  from  this  increment  of  capability  is 
0.8  x  1 3  =  10.4,  or  10.4  =  ki^2i  units.  On  the  third  increment 
of  30  units  (22,001  to  23,000  seats)  the  expected  revenue  is 
0.73  x  13  -  3.8  =  ki<f31  units  per  unit  of  capability  since  there 
is  a  23  per  cent  chance  that  none  of  these  units  of  capability 
will  be  used  and  73  per  cent  that  all  will  be  used.  For  the 
fourth  Increment  of  20  units  (23,001  to  27,000  seats)  of  capability 
the  expected  revenue  Is  0.4  x  13  =  3-2  =  kitf41  units  per  unit  of 
capability,  while  for  the  fifth  Increment  of  30  units  (27,001  to 
30,000  seats)  it  is  0.2  x  1 3  =  26  =  k^i  units  per  unit.  For 
the  sixth  increment,  which  is  the  number  of  units  assigned  above 
the  30,000  seat  mark,  the  expected  revenue  is  0.0  x  13  =  0  per 
unit  since  it  is  certain  that  none  of  these  units  of  capability 
can  be  used.  It  is  clear  that  no  assignments  above  30,000  seats 
are  worthwhile,  and  hence  the  last  increment  can  be  omitted. 
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The  index  h  =  1,  2,  ),  4,  5  will  be  used  to  denote  the  l3t,  2nd, 

. ..,  ljth  increment  of  demand. 

The  number  of  assigned  units  in  each  increment,  however, 
can  be  viewed  as  an  unknown  that  depends  on  the  total  (passenger- 
hauling)  capability  assigned  to  Route  J  -  1 .  Thus  if  the  total 
assigned  is  Yi  =  210  units  of  capability  then  the  part  of  thi3 
total  belonging  to  the  first  increment,  denoted  by  yM,  is 
y i i  -  200  and  the  part  belonging  to  the  second  Increment,  de¬ 
noted  by  y2  j  ,  is  y21  -  10;  tiie  amounts  in  the  higher  increments 
are  y^  -  0  for  i  -  j>,  4 ,  g .  To  review,  the  passenger-carrying 
capability  Y,  is  determined  by  the  number'  of  aircraft  assigned 
to  Route  ,j ,  so  that 


b) 


Yj  -  pu  xu  +  p 


2J  X2J  +  P^J 


“JJ 


+  Ph 


V4J 


On  the  other  hand,  Y,  itself  breaks  down  into  five  increments 


(°>  y|  '  yU  '  y2j  +  gj  "  +  g.l 


for  Routes  J  =  1,  J,  4,  and  correspondingly  fewer  for  J  -2,  5* 


Regardless 

of  the  total 

V 

the  amount 

J'hj 

belonging  to  each 

increment  ' 

Ls  bounded  by 

the 

total  size 

b,  , 
dj 

of  that  increment; 

the  latter,  however,  is  simply  the  change  in  demand  level,  so 
that 


0 

< 

yu 

< 

d,  . 

1 J 

-  V 

0 

/ 

y2j 

< 

d2j  ' 

_  dlJ  : 

=  b2  j  ’ 

0 

< 

/ 

"  d2J  ' 

0 

< 

y,u 

< 

d4j  - 

-  d  ,  .  - 

'  b4  J  ’ 

0 

/ 

y..  . 

< 

dr  .  - 

-  d,  .  = 

-  b  .  . 

— 

j  J 

4  J 

J  J 

(7) 
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The  total  expected  revenue  from  Route  J  is,  therefore, 

(8)  kj(^J  ylJ  +  ^2J  y2J  +  +  *9J  y5J  ^ ' 

where  k,  is  revenue  (in  thousands  of  dollars)  per  100  passengers 
J 


carried 

on  Route  J, 

and 

where 

,  as  s 

een  in 

Table 

(9) 

1  *  bj  - 

hj 

+  X2J 

+  Lj 

+ 

+ 

^  - 

A2J 

+ 

+  *i,  j 

+ 

^5J * 

li 

>^» 

Aa 

+ 

4- 

Cij  - 

A,u 

+ 

Lj  - 

AbJ‘ 

For  example,  the  total  expected  revenue  for  Route  1  is 

(10)  1  ( 1  *0y i !  +  .8y 1 2  +  -Thy  13  +  %  1 4  +  •  2yis). 


The  most  important  fact  to  note  about  the  linear  form  (10)  i3 
the  decrease  in  the  successive  values  of  the  coefficients  y,  . . 

hj 

Moreover,  this  will  always  be  the  case  whatever  the  distribution 
of  demand  since  tne  probability  of  equaling  or  exceeding  a  given 
demand  level  d,  .  decreases  with  increasing  values  of  demand. 

Suppose  now  that  ytl,  y2i,  ...  are  treated  as  unknown 
variables  in  a  linear-programming  problem  subject  only  to  (6) 
and  (7),  where  the  objective  is  to  maximize  revenues.  Let  us 
suppose  further  that  Yi  is  fixed.  It  is  clear,  since  the  co¬ 
efficient  of  y 1 !  is  largest  in  the  maximizing  form  (8),  that 
ytl  will  be  chosen  as  large  as  possible  consistent  with  (6)  and 
(7) ;  for  the  chosen  value  yn,  the  next  increment  y2 1  will  be 
chosen  as  large  as  possible  consistent  with  (6)  and  (7),  etc. 


hm-1  f)y$ 
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Thus,  we  need  only  specify  y^  by  restrictions  (6)  and  (7), 
because  when  the  maximum  Is  reached  the  values  of  the  variables 
y  m  ,y2 1  ,  ...  are  precisely  the  Incremental  values  (6)  associated 
with  Yj.  Even  if  passenger  capability  Yj  is  not  fixed,  as  in 
the  case  about  to  be  considered,  it  should  be  noted  that  what¬ 
ever*  be  the  value  of  Yj  the  values  of  yn,y2i  ...  that  minimize 
an  over-al.l  cost  form  such  as  (l4)  below  must  maximize  (8)  for 
,i  -  1  ,  so  that  the  incremental  values  of  Y,  will  be  generated 

by  y 1 1 #  y  2 1  >  •  •  •  • 

The  1  i ne a r-p rot; ramming  problem  In  the  case  of  uncertain 
demand  becomes  : 

UNCERTAIN  DEMAND  MODEL 

Find  numbers  x,  .  and  y,  . ,  and  the  minimum  value  of  z. 
ij  LJ 

satisfying  the  following  conditions. 

(11) 

(12) 

(U) 

(14) 


Row  Sums  : 

Xil 

+  X 

12 

•f 

•  • 

+  xin 

=  ai 

(i  = 

i, 

2 ,  .  .  . , 

m) 

Column 
Sums  : 

PU- 

X1J 

+  P 

2JX 

2J 

\  ... 

H  pmj 

X  1 

nj 

= 

yu 

+ 

y2J 

+ 

. .  .  + 

yrj 

(J  =  1, 

2, 

•  ■ « ,  n- 

-1) 

*lj 

>  0 

f 

(i 

=  1, 

.  .  . ,  m 

J  J  =  1 

) 

.  n) 

o 

l/\ 

yhJ 

< 

b.  , 
hj 

(h  =  1; 

i  •  *  *  i 

r;  J  - 

1, 

•  •  • ,  n- 

-i) 

Expected 

m 

n 

- 

n-1 

r 

- 

Costs  : 

Z  — 

i 

1=1 

l 

J  = 

c 

1 

ij 

X1J  + 

Ro- 

E  ki 

j=i  J 

I 

h=l 

yhJyhJ 
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lie  re  Rq  In  the  value  that  expected  revenue  would  be  If  sufficient 
aeat3  were  supplied  for  all  customers.  Thu3  expected  costs  are 
defined  as  total  outlays  (first  term)  plus  the  expected  loss  of 
revenue  due  to  shortage  of  seats  (last  two  terms). 

For  the  problem  at  hand,  the  bounds  b^  and  the  expected 
revenues  per  unit  for  the  "incremental  variables"  y^j  can 
be  computed  from  the  probability  distributions  of  Table  6  via 
(7)  and  (9) . 

The  numerical  values  of  the  constants  for  the  stochastic 
case  are  shown  in  Table  7* 


Table  7 


INCREMENTAL  ROUNDS  bhJ  AND  EXPECTED  REVENUES  kj 
PER  UNIT  OF  ASSIGNED  PASSENGER-CARRYING  CAPABILITY 


Route  1 

Route  2  j 

Route  5 

Route  4 

Route  5 

Incre¬ 
ment  h 

b, 

hi 

bM 

bh5 

Rr 

5  h5 

1 

200 

15 

50 

15 

140 

7 

10 

7 

580 

1 

O 

20 

BSD 

100 

9.1 

20 

6.5 

40 

5  -6 

20 

0.9 

5 

50 

9.8 

*  * 

20 

4.9 

30 

4  .2 

20 

0.1 

A 

20 

5.2 

X  -If 

20 

2.1 

20 

.  1 

2.1 

** 

5 

30 

2.0 

*  -H 

2C 

0.7 

240 

0.7 

*# 


Only  two  increments  for  Route  2  and  three  increments  for 
Route  5  are  needed  to  describe  the  distribution  of  demand. 
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4 .  RULES  FOR  COMPUTATION 

The  work  sheet  for  determining  the  optimal  assignment 

under  uncertain  demand  is  shown  in  Table  9.  To  form  the  new 

row  equations  (ll),  the  x,  ,  entries  are  summed  to  yield  the  a. 

1  J  1 

values  given  in  the  Aircraft  Available  column.  To  form  the 
column  equations  (12),  the  ,  entries  are  multiplied  by  the 
corresponding  p.  ,,  the  y,  ,  by  -1,  ana  summed  down  to  yield  zero. 

Step  1 .  To  initiate  the  computation  any  set  of  non-negative 
values  may  be  assigned  to  the  unknowns  x^  ,  and  y^  provided  they 
satisfy  the  equations  and  thereby  constitute  a  feasible  solution. 

Step  2  .  Circle  any  m  +  n  of  the  x^  ,  and  y^  entries,  where 
m  +  n  is  the  total  number  of  row  and  column  equations.  These 
circles  can  be  arbitrarily  selected  except  that  they  must  have 
the  property  that  if  the  fixed  values  assigned  to  the  other  non- 
clrcled  variables  and  the  constant  terms  were  arbitrarily  changed 
to  other  values  then  the  circled  variables  would  be  determined 
uniquely  in  terms  of  the  latter.  Such  a  circled  set  of  variables 
is  called  a  basic  set  of  variables:  the  array  of  coefficients 
associated  with  this  set  in  the  equations  (l))  and  (12)  is  re¬ 
ferred  to  as  the  basis  in  the  theory  of  the  simplex  method  [4]  . 
Note  :  One  simple  way  of  selecting  a  basic  set  is  shown  in 

Table  10.  One  x.  .  entry  is  arbitrarily  selected  and  circled  in 

1J 

each  row  corresponding  to  a  row  equation,  and  one  y^,  is  arbi¬ 
trarily  selected  and  circled  in  each  column  corresponding  to  a 
column  equation  In  general,  it  is  suggested  that  entries  be 
circled  that  appear  to  have  a  chance  of  having  a  positive  value 
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in  an  optimum  solution;  for  values,  the  last  entry  in  the 
column  that  appears  likely  to  be  positive  in  an  optimum  solution 
should  be  circled. 

Step  j.  For  (i,J)  and  (h,j)  combinations  corresponding  to 
circled  entries,  compute  implicit  prices  u^  and  v^  associated 
with  equations  by  determining  values  of  u^  and  Vj  satisfying  the 
equations 


(15) 

ui  +  piJvJ  0  ciJ 

(x^j  circled), 

(16) 

0  +  (-l)Vj  =  -kj»hJ 

(yhJ  circled)  . 

There  are  always  m  +  n  equations  (15)  and  (16)  in  m  +  n  unknowns 
u.^  and  v^  that  can  be  shown  easily  to  have  a  unique  solution  [4]  . 
They  can  be  solved  by  inspection,  for  it  can  be  shown  that  the 
system  either  is  completely  triangular  or  at  worst  contains  sub¬ 
systems — some  triangular  and  some  triangular  if  one  unknown  is 
specified . * 

Step  4 .  For  each  box  corresponding  to  x^j  or  y^^,  compute 


(17) 

6lj  -  <ui  +  pljvj>  -  C1J 

(for  x^j  box) , 

(18) 

6hJ  '  <°-  vj>  -  '-k/hj) 

(for  yhJ  box)  . 

*Thls  is  the  analogue — for  the  "generalized"  transportation 
problem  (1),  (2),  (5),  (4) — of  the  well-known  theorem  for  the 
standard  transportation  problem  that  all  bases  are  triangular. 
Its  proof  is  similar. 
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In  practice,  one  of  the  6,  .  or  8'  is  recorded:  the  others  are 

ij  hj 

computed  and  compared  with  it,  and  the  largest  in  absolute  value 

is  used.  It  can  be  shown  Pil  that  if  the  x.  ,  or  y,  .  value  as- 

L  J  ij  Jhj 

sociated  with  a  noncircled  entry  is  changed  to 

+  0  or  yh  !  ±  6  (©  >  0) , 

the  other  noncircled  variables  remaining  invariant,  and  the 
circled  variables  adjusted,  then  the  expected  costs  z  will  change 
to  z 1 ,  where 


z ' 

cr 

z  +  96^  or 

z '  =  z 

+  eSg  . 

Thus 

it  [jays 

to 

increase  x^j  or 

yhj 

if  6± j  or 

^hj  y  °’ 

unless 

yhj 

is  equal 

to 

its  upper  bound 

bh.J  * 

in  which 

case  no 

increase 

in  is  allowed;  also  it  pays  to  decrease  x^ (  or  y^j  if  8^ 
or1  8^j  <  0  unless  x^  (  =0  or1  y^j  =  0,  in  which  case  no  decrease 
is  allowed  , 

Test  for  Optimality:  According  to  the  theory  of  the  simplex 
method  [y]  if  the  noncircled  variables  satisfy  the  following 
conditions! 

(a)  each  one  is  at  either  its  upper  or  its  lower  bound 
value , 

(b)  the  corresponding  5^.  or  b^  ,  is  less  than  or  equal 
to  0,  if  it  is  at  its  lower  bound  value,  and 

(c)  the  corresponding  8^  ^  or  5^  .  is  greater  than  or  equal 
to  0  if  it  is  at  its  upper  bound  value, 

then  the  solution  is  optimal  and  the  algorithm  terminates  .  Other¬ 
wise  there  are  5^  ,  or  5^  .  for  which  a  decrease  or  increase 
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(depcndlng  on  whether  the  sign  is  negative  or  positive)  in  the 

corresponding  variable  is  allowed;  let  the  largest  among  them 

in  absolute  value  be  denoted  by  &  or  8 1  . 

J  rs  rs 

Step  5-  Leaving  all  nonclrcled  entries  fixed  except  for 
the  value  of  the  variable  corresponding  to  the  (r,s)  box  deter¬ 
mined  in  Step  '( ,  modify  the  value  of  x  (or  y  )  to 

rs  ps 

xr3  +  Q  (or  yr3  +  0)  if  Srs  >  0  (or  >  0) 

xrs  ~  0  (or  yrs  ~  if  6rs  <  0  (or  6rs  <  °)' 

where  0  >  0  is  unknown,  and  recompute  the  values  of  circled 

variables  as  linear  functions  of  0.  Choose  the  value  of  ©  =  9* 
a3  the  largest  value  possible  consistent  with  keeping  all  basic 
(circled)  variables  [whose  values  now  depend  on  0]  between  their 
upper  and  lower  bounds;  in  the  next  cycle  correct  the  values  of 
the  circled  variables  on  the  assumption  that  0  =  0*. 

Also,  if  at  the  value  0=0*  one  (or  more)  of  the  circled 
variables  attains  its  upper  or  lower  bound,  in  the  next  cycle  drop 
just  one  such  variable  from  the  basic  set  and  circle  the  vari¬ 
able  x  instead.  Should  it  happen  that  it  is  x  that  attains 
rs  rs 

its  upper  or  lower  bound  at  0  =  0*,  the  set  of  circled  variables 

is  the  same  as  before;  their  values,  however,  are  changed  to 

allow  x  to  be  fixed  at  its  new  bound, 
rs 

Start  the  next  cycle  of  the  iterative  procedure  by  return¬ 


ing  to  Step  5  • 
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S  NUMERICAL  SOLUTION  OF  THE  ROUTING  PROBLEM 

For  our  storting  solution  we  used  for'  values  of  the  x^ 
the  best  solution  ol’  the  earlier  study,  assuming  fixed  demands 
equal  to  the  expected  value’s  of  the  distribution.  These  are 
shown  in  Table  10.  These  x^  ,  will  meet  the  expected  demands, 
so  that  Y.  b.  except  for  Route  U;  there  is  a  deficit  of  100 
1'or’  this  route,  and  by  (b)  we  have  Y3  -  '  00.  These  Y,  are  broken 

J 

down  into  the  successive  incremental  values  sliown  below  the 
double  line  in  Table  10;  see  Eq .  (6). 

Next,  one  of  the  variables  in  each  row  is  circled.  In  the 
example,  the  selected  variable'  are  x^,  x22,  x35,  and  x43;  each 
appears  likely  to  be  in  an  optimal  solution,  though  x4 x  may 
turn  out  to  be  a  better  choice  than  x43.  Next,  the  last  positive 
entry  in  each  column  is  circled;  in  the  example,  these  are  the 
variables  y31,  y22,  ya3,  y-w,  and  y15.  In  all,  there  are  m  +  n 
circled  variables  (9  in  the  example).  The  implicit  values  must 
satisfy  the  m  +  n,  or  9,  equations: 


*  II 

In  the  humorous  parody  by  Paul  Gunther,  entitled  Use  of 
Linear  Programming  in  Capital  Budgeting,"  Journal  of  the  Opera¬ 
tions  Research  Society  of  America,  May,  1935;  Mrs.  Efficiency 
wondered  why  Mr.  0.  R.  did  not  start  out  with  a  good  gue3s.  It 
will  be  noted  that  in  this  paper  we  have  followed  Mrs.  Efficiency's 
suggestion  and  have  started  with  a  guess  at  the  final  solution 
rather  than  going  through  the  customary  use  of  artificial  vari¬ 
ables  and  a  Phase  1  of  the  simplex  process. 
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Ui 

p  1  1  v  1 

c  1  1 

D 

11 

a. 

c  1 1 

-18), 

u2 

-t  p22V2 

-  C22 

(p22=10, 

c22 

=  18), 

u  3 

+  P  3  5  V  ‘5 

=  C  3  5 

(P35=29 , 

C35 

=6), 

U4 

+  p43V3 

=  c4  ^ 

(p43=22 , 

C43 

=  17), 

0  f 

(-1  )v. 

=  — k  1  ^3  1 

(k  1=9 • 

8), 

0  + 

(-1 ) V2 

=  —  k2^22 

(k2^22=9 • 

1), 

0  + 

(-1 )  v3 

=  — k3^33 

(k3*33=4. 

9), 

0  + 

(  1  )  V  4 

=  —  k  4  ^4  4 

( k4^ 44=2 . 

1), 

0  + 

(-1 )vs 

=  -i<  5  y  1 5 

( k5^ 15=1. 

0). 

permits  the 

computation  of  8^ 

and  8  ’  . ;  3ee 

(17)  and 

As  a  check,  8^  =  0  and  8^  =  0  for  (i,J)  and  (h,J)  corresponding 
to  circled  variables.  The  6^  or  of  largest  absolute  value 
is 


b24  =  [-76  +  16  ( 2 . 1 )]  -1 4  =  -08  .6; 


hence  a  decrease  in  the  variable  x24  with  adjustments  of  the 
circled  variables  will  result  in  a  decrease  in  the  expected 
costs  by  an  amount  of  58-5  units  per  unit  decrease  in  x24.  If 
x24  =  b  is  changed  to  x24  =6-9,  then,  in  order  to  satisfy 
the  column  4  equation,  the  circled  variable  y44  =  10  must  be 
modified  to  y44  =  10  -  160  (all  other  variables  in  column  4  are 
fixed).  Also,  to  satisfy  the  row  2  equation,  x22  *  8  must  be 
modified  to  x22  =8+9;  this  in  turn  causes  y22  =  70  to  be 
changed  to  y22  =  70  +  109  in  order  to  satisfy  the  column  2 
equation.  The  largest  value  of  9  is  9*  =  1 0/ 13/  at  which  value 

y44  =  0. 
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The  numerical  valued  o!'  the  variables  appearing  In  Table  11 
ir'.‘  obbalnei  f i’o.t.  those  of  Table  10  by  setting  Q  -  0*  -  10/lS. 

Trie  v  rlab  l*  *  x24  becomes  :  new  :  Ire  led  variable  in  place  of  y44 , 
which  nil  Ilr.  lower  bound,  " ' ' r ■  :  the  other  variables  to  be 
circled  remain  the  same  as  In  Table  10.  Computing  the  new  set 
of  lrr;  11  'll  prices,  we  see  that  the  5^  t  of  largest  absolute  value 
tnat  wan  bo'  Increased  or’  decreased  (accor'diny,  to  the  siryn  of  5,.) 
is  S2 3  2MM  Chan.-.ln;  x23  to  -  0  requires  that  the  variables 
*22  >  y a;.-,  and  Yss  be  modified  is  shown  in  Table  11.  The  maximum 
value  of  0  Is  0  0"  -  20/M,  at  which  value  we  have  y33  =  0. 

The  new  solution,  In  which  x23  replaces  y33  os  a  circled  variable, 
l.s  ,r,l  ven  In  Tabl“  12,  whor>*  the  dccr’ease  in  the  noncircled  vari¬ 
able  x,t  j  causer,  chamyes  In  the  variables  x43,  x22  ,  x23 ,  y3], 
and  y22 .  The  Jar.yest  value  of  0  Is  j/lo,  at  which  value  y22 


hits  Its  upp-.’r  bound  b22  100. 

In  the  pa.-.so.'/e  from  Table  1)  to  Table  M  we  have  taken  a 


"double" 

s  t  e  p  . 

The  maximun.  Increase 

Is 

0 

-  8o/20,  at  which 

point 

.Via  alts 

ItS  U| 

;  e  r  hound  b \ 3  -  ■  0  . 

1 1 

Is 

easy  to  see  that 

if 

ri  :xt  the 

Lnr  !’•  •; 

i.ental  variable  y23  is 

In  : 

re 

nse.i  tnen  8 ass 

oc  i  at*: 

w  i  t  n  x  3  ;j 

snou  1 

1  1)  '  SI  lllj"  i  to  8  •)  ;»  + 

1  ( 

a  ^2 a  ) *< 5  -  —'My 

4- 

2  •(  .O  -  •  ’)  .a;  ther-fot*.' ,  It.  Is  economical  to  Increase  y2i> 

is  wall  as  y  j .  ilowo’ver’,  it  can  b  '  shown  that  the  sl;yn  of  S32 
woulu  become  positive  if  the  next  increment,  y3^,  were  considered. 
The  maximum  value  of  0-0*  is  100/2. >. 

In  the  passage  from  Table  12  to  ,  It  will  be  noted  that 
the  variable*  y33 ,  wni.an  na  i  been  arop.ped  earlier,  is  a.-yuln  b roubinl 


Into  the  solution.  The  maximum  value  of  0  is  22/20,  at  which 
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value  y33  reaches  its  upper  bound,  so  that  the  new  solution, 
given  in  Table  13,  has  the  same  set  of  circled  variables  and 
hence  the  same  implicit  values  as  those  in  Table  14 .  Moreover, 
the  solution  is  optimal  since  all  noncircled  variables  are  either 
at  their  upper  or  lower  bounds — those  at  upper  bounds  have  cor¬ 
responding  5^  ,  >  0  and  those  at  lower  bounds  have  6^ j  <0* 

In  comparing  thl3  solution  (Table  19)  with  the  optimal 
solution  for  the  fixed-demand  case  (Table  3),  it  is  interesting 
to  note  that  the  chief  difference  appears  to  be  a  general  ten¬ 
dency,  in  the  case  of  distribution  with  sharp  peaks,  to  shift 
the  total  seats  made  available  on  route  to  a  mode  of  the  dis¬ 
tribution  rather  than  to  the  mean  of  the  distribution.  The 
total  3eats  made  available  on  routes  with  flat-  distributions 
of  demand,  on  the  other  hand,  appear  to  be  at  the  highest  level 
attainable  with  the  residual  passenger-carrying  potential. 

To  compute  the  expected  cost3  of  the  various  solutions, 
the  first  step  (see  Eq .  (14))  is  to  determine  what  the  expected 
revenues  Rq  would  be  if  sufficient  seating  capacity  were  furnished 
at  all  times  to  supply  all  passengers  that  show.  From  Table  2 
it  Is  easy  to  see  tnat 

RQ  =  13(230)  +  13(120)  +  7 ( 1 80 )  +  7(90)  +  1(600)  =  7300, 
so  that  the  expected  revenue  would  be  $7,300,000. 
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Table  8 


COMPARATIVE  COSTS  OF  VARIOUS  SOLUTIONS 


Table 

Expected 

Revenues 

For 

Seats 

Supplied 

(1) 

Expected 

L03t 

Revenues* 

(2) 

Operating 

Costs 

(3) 

Net 

Expected 

Cost 

(Thousands ) 

(2)  +  (8) 

10 

-odd1* 

766 

900 

1 , 666 

1 1 

-6674 

72o 

901 

1,627 

12 

-<1607 

698 

901 

1,694 

'•8 

—6688 

o62 

899 

1  ,pbl 

14 

— o64  l 

669 

888 

1,64  2 

1  'J 

— oG^  9 

o4  i 

0  0 
uu) 

1,642 

Data  In  column  (2)  are  obtained  by  subtracting  the 
expected  revenues  for  seats  supplied,  column  (l),  from  R»  =  7800 
~  tne  expecteu  revenues  II’  an  unlimited  number  of  seats  were 
supplied  . 


it  is  seen  tn  it  the  solution  presented  in  the  earlier  paper 
[l]  ,  assuming  demands  to  be  exactly  equal  to  their  expected  values, 
has  a  net  expecteu  cost  of  $’.,uoO,000.  [it  is  interesting  to 
note  that  if  the  demands  were  fixed  and  equal  to  their  expected 
values,  the  costs  would  be  only  $1,000,000  (see  Table  8).  The 
o7  per  cent  increase  in  net  cost  for  the  variable-demand  case 
is  due  to  18,400  additional  passengers  (on  tne  average)  being 
turned  away  because  of  the  distributions  of  demand  assumed  in 
Table  6.]  The  successive  improvements  in  the  solution  given  In 
Tables  10  to  1;  reduce  the  net  expected  costs  from  $1,666,000 
to  $1,  2h,000  for  tne  optimal  solution. 
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In  the  illustration  the  best  solution  obtained  by  pretend¬ 
ing  that  lemands  are  fixe  i  at  th eye  expected  values  has  a  9  per 
cent  higher  expected  rest  than  that  for  the  best  solution  obtained 
by  usimj;  the  assumed  distributions  o!'  demand.  It  is  also  seen 
that  very  little  additional  reputational  effort  was  required 
to  take  account  of  this  uncertainty  of  demand. 


Table  9 
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.v'CRK  SHEET  FOR  DETERMINING  OPTIMAL  ASSIGNMENT  UNDER  UNCERTAIN  DEMAND 


Bex  not  used  because  correspond inr;  row  or  column  has  no  equation, 
or  else  because  airernPt  type  cannot  Ply  required  ran^e ,  or  fewer  in¬ 
crements  are  needed  to  describe  the  distribution  of  demand  on  the  route. 


Table  10  -  Cycle  0 
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WORK  SHEET  FOR  DETERMINING  OPTIMAL  ASSIGNMENT  UNDER  UNCERTAIN  DEMAND 
624  =  3  8 .  A  ,  9  -  lO/lL,  Expected  Coat  =  $1  ,060, 000 


Route 


0) 

N.  Y.  N.  Y.  N .  Y 

to  to  to 

L.  A.  L.  A.  Dali 

1-3  top  2-stoo  0-st 


(1)  A 


(2)  B 


(O  C 


CO  D 


Incre-  200 
mcnt 

(,)  -e 


CO  (6) 

N.  Y.  Sur-  Air- 

to  plus  craf 

Boston  Air-  Aval 

0-stop  craft  able 

8l  0 


-1 
-10. A 


s 

|  3° 

-1 

-1 

-A  .9 

-A  .2 

0-ise 

-l 

-l 

-2.1 

-2.1 

Air¬ 

craft 

Avail¬ 

able 

Im¬ 
plicit 
Prices 
.  ui 

10 

-139 

19 

1 

-  7o 

Net 


Im¬ 

plicit 

Prices 


v  , 


-1 

-2  .6 


0 
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_^4_  Table  11  -  Cycle  1 

WORK  SHEET  FOR  DETERMINING  OPTIMAL  ASSIGNMENT  UNDER  UNCERTAIN  DEMAND 
=  2}.4,  0  =  20/14,  Expected  Cost  =  $1,627,000 
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Table  12  Cycle  2  -35- 

20 HK  SHEET  FOH  DETERMINING  OPTIMAL  ASSIGNMENT  UNDER  UNCERTAIN  DEMAND 
641  =  -50.8,  0  -  9/1 -j.  Expected  Cost  =  $1,59*1,000 


Route 


0)  (2)  a)  ('•)  (b) 

Type  N.  Y.  N.  Y.  N.  Y.  N.  Y.  N.  Y. 

of  to  to  to  to  to 

Air1—  L.  A  L.  A.  Dellas  Dallas  Boston 

craft  ]— stop  2— stop  0-stop  1-stop  0— stop 


able 


10 


Im¬ 

plicit 

Prices 

u. 


HM-1  cro 
1 2-(-  b 

_j()_  Table  1  j  -  C; 

.•10 RK  .  MEET  FOR  DETERMINING  OPTIMAL  AS: 
S32  =  b  ■ :  ,  ©  =  1 00/29  -  ,  E. 


Route 

(:) 

(*) 

(  5) 

('• 

Type 

M.  Y. 

N.  V. 

M.  Y. 

N.  Y 

o  r 

to 

to 

to 

to 

A  1  r- 

Ij  .  A  . 

L .  A  . 

Da  J  1  a  .a 

!  Dali; 
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'I  1  i >  1  f '  H  -  Cycle  ;i  —57— 

WORK  SHEET  FOR  DETERMINING  OPTIMAL  ASSIGNMENT  UNDER  UNCERTAIN  DEMAND 
fi-j 3  :  -.f»,  0  -  20/22  -  .9,  Expected  Cost  =  $1,542,000 
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_j,3_  Table  15  -  Cycle  5 

(Optimal) 

WORK  SHEET  FOR  DETERMINING  OPTIMAL  ASSIGNMENT  UNDER  UNCERTAIN  DEMAND 
_  Minimum  Expected  Cost  11,524,000 

Route 


(1)  (2)  (5)  (4 

Type  N.Y.  N.Y.  N.Y.  N. 

of  to  to  to  •  to 

Air-  j L .  A.  L.  A.  Dallas  Dal 

craft  I 1-stop  2-stop | 0-stop  1-s 


Air-  Im- 

craft  plicit 
Avail-  Prices 
able  u 


1 

10 

0  1 

0  ! 

km 

1 3-7-yo 

-yy 
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